Disorder Induced Phases in Higher Spin Antiferromagnetic Heisenberg Chains 



Enrico Carlon/'Q Peter Lajko, 2,3 Heiko Rieger, 1 and Ferenc Igloi 3,4 

1 Theoretische Physik, Universitat des Saarlandes, D-66041 Saarbrucken, Germany 
2 Department of Physics, Kuwait University, P. O. Box 5969, Safat, Kuwait 13060 
' Institute of Theoretical Physics, Szeged University, H-6720 Szeged, Hungary 
^Research Institute for Solid State Physics and Optics, H-1525 Budapest, P.O.Box 49, Hungary 

(Dated: February 2, 2008) 

Extensive DMRG calculations for spin S = 1/2 and S — 3/2 disordered antiferromagnetic Heisen- 
berg chains show a rather distinct behavior in the two cases. While at sufficiently strong disorder 
both systems are in a random singlet phase, we show that weak disorder is an irrelevant perturbation 
for the S = 3/2 chain, contrary to what expected from a naive application of the Harris criterion. 
The observed irrelevance is attributed to the presence of a new correlation length due to enhanced 
end-to-end correlations. This phenomenon is expected to occur for all half-integer S > 1/2 chains. 
A possible phase diagram of the chain for generic S is also discussed. 

PACS numbers: 05.50.+q, 64.60.Ak, 68.35.Rh 



I. INTRODUCTION 

A simple model where the effects of the interplay be- 
tween quantum fluctuations and disorder can be studied 
in great detail is the spin-5 random antiferromagnetic 
Heisenberg (AFH) chain with Hamiltonian: 

i 

where the Ji > are quenched random variables. 

In absence of randomness (Ji = J) quantum fluctu- 
ations lead to qualitatively different behavior for half- 
integer (5 = 1/2, 3/2 . . .) and integer (S = 1,2...) val- 
ues of the spin Half- integer spin chains have a gap- 
less spectrum and quasi long range order. It is believed 
that they all belong to the same (bulk) universality class 
independently on 5 01 ■ This was explicitly verified nu- 
merically for the 5 = 3/2 chain which was found to 
have the same bulk decay exponent as for the 5=1/2 
chain. Integer spin chains are instead gapped and have 
a hidden topological order 

The effect of quenched randomness on the AFH chains 
can be studied by the Strong Disorder Renormalization 
Group (SDRG) method. This technique consists in suc- 
cessively decimating out the strongest bond in the chain 
H. For the 5 = 1/2 chain the SDRG procedure can be 
carried out in great detail and yields a probability distri- 
bution of coupling constants that under renormalization 
is broadened without limits and in the so called infinite 
randomness fixed point the method becomes asymptoti- 
cally exact. At this fixed point, describing the so-called 
random singlet (RS) phase [||, the energy scale (r) and 
the length scale (£) are related via lnT ~ which dif- 
fers from scaling at a conventional (random) fixed point, 
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where T ~ £~ z with a finite dynamical exponent z. Nu- 
merical studies on random spin 5 = 1/2 AFM chains 
generally agree with the SDRG results 0, , although 
some issues are still debated [jj. 

For 5 > 1/2 the SDRG procedure requires a higher de- 
gree of approximations with respect to the 5=1/2 case. 
A spin 5 is represented by a maximally symmetrized 
combination of 25 identical 5=1/2 spins. The SDRG 
decimation then produces effective spins of magnitudes 
5 c ff < 5 obtained by linking the spin-1/2 objects into 
correlated singlets. For the 5=1 case an SDRG analysis 
shows that a sufficiently strong randomness induces a RS 
phase with S e ff = 1, separated from the gapped phase by 
a gapless region of Griffiths singularities |lfj (the gapless 
Haldane phase 11]) in which the dynamical exponent, 
z, increases with disorder. A recent SDRG study 0] °f 
the disordered AFH chain with 5 = 3/2 indicated the 
existence of two phases: At strong disorder the relevant 
degrees of freedom are 5 e ff = 3/2 spins, while at weaker 
disorder they are of type S e g = 1/2. Both phases are 
of RS type with identical critical exponents. A quantum 
critical point separating them is expected to have specific 
multicritical exponents [T2I H"^ . 

As the SDRG method might fail at weak disorder, its 
predictions should be tested by means of accurate numer- 
ical analysis, which is the aim of the present paper. We 
will present numerical evidences, supported by theoreti- 
cal arguments, that the phase diagram for the 5 = 3/2 
spin chain differs from what predicted by SDRG [l2fl . 
Our numerical results support the existence of a RS phase 
only for sufficiently strong disorder, while weak disorder 
appear to be an irrelevant perturbation for the system. 



II. NUMERICAL RESULTS 

Here we analyze the critical behavior of the AFH 
chains with varying degree of disorder by means of den- 
sity matrix renormalization group (DMRG) techniques 
jl4J ] . Numerical calculations are restricted to the 5 < 3/2 




FIG. 1: Average correlations in the random S = 1/2 and 3/2 
chains with disorder strengths 8 — 0.3 (circles) and S = 0.8 
(squares), (a) End-to-end correlations; dashed lines are linear 
fits of the data.(b) Cl vs L in a log-log scale. Here the dashed 
line show the decay in the pure system (1/L) and in a RS 
phase 1/L 2 . 



FIG. 2: Integrated probability distribution of the first gap in 
log-log scale, (a) S = 3/2, 5 = 0.3 (L = 4,8,12,16,20,24,32 
from left to right) and 5 = 3 (L = 4,6,8,10) shown as solid 
lines, (b) 5" = 1/2, S = 0.3 (L = 4,. . . ,24,32) shown as dashed 
lines. Inset: Collapse of the data for S = 3/2 and S — 3 as 
expected in a RS phase (L — 4, 6, 8, 10). 



case, and the disorder average is performed on 10 3 - 10 4 
samples taken from the distribution: 

PS (J) = 5- 1 J- l + l / s for < J < 1 (2) 

where <5 2 = var[ln J\ measures the strength of disorder. 
As usual in DMRG LJ we use open boundary conditions 
and keep typically m — 80 — 100 states in the renormal- 
ization procedure. In order to obtain numerically stable 
results we use the DMRG finite system method perform- 
ing several "sweeps" (~ 4 — 5) through the chain [l4j . 
which is an essential step in the DMRG procedure. The 
DMRG calculations for the S = 3/2 chain have been ex- 
tended up to L = 32 for weak disorder, while at stronger 
disorder the procedure tends to become less stable and 
we had to restrict ourselves to shorter chains. However 
the physical behavior in the strong disorder regime turns 
out to be quite well understood, while, in view of the 
limitations of the SDRG at weak disorder, it is precisely 
this regime which is physically the most interesting. In 
the weak disorder regime the DMRG is a rather stable 
and reliable technique. 

A. Correlation functions 

In order to emphasize the differences between the two 
cases we present together the results for the 5 = 3/2 
and S = 1/2 chains for the spin-spin correlation function 
defined as: 

C(i,j) = [(S!S?)U (3) 

where [. . .] av denotes the averaging over quenched disor- 
der. The correlation function obviously does not depend 
on the spin direction. 



Figure [Ha) shows a plot of the surface-surface corre- 
lation function C£ = C(1,L) plotted as a function of 
1/L for S =1/2(5 = 0.3) and S = 3/2 (5 = 0.3, 0.8). 
In a RS phase one expects that asymptotically in L the 
surface-surface correlation vanishes as C S L ~ L~ v ks with 
Vrs = 1 8]. This is indeed observed in the case S = 1/2. 
In the S — 3/2 case instead C£ extrapolates to a finite 
value which indicates the presence of a surface ordering 
phenomenon or, in other words, a first order surface tran- 
sition (formally r) s =0). The same type of behavior is 
observed for the 5' = 3/2 chain in absence of randomness 


To obtain the bulk exponent (77) we consider C(l, L/2) 
the correlation function between an edge and a spin in 
the middle of the chain, which decays asymptotically as 
C(l,L/2 + 1) - L~(v+v')/z, Thus each point of refer- 
ence brings its local scaling dimension, r\/2 and t] s /2, 
respectively. In order to eliminate the surface exponent 
contribution we consider the ratio 

_ C(l,L/2 + l)C(L/2,L) 

= cl ^ ' 

which decays asymptotically as Cl ~ L~ v . The plot of 
lnC^ vs. InL is shown in Fig. ^b). In the pure system 
the bulk exponent is r\ = 1, while in an RS phase one 
expects tjyis — 2. The numerical data for S = 3/2 are 
quite consistent with a 1/L decay, which seem to be in 
disagreement with the existence of a RS phase at weak 
disorder, but rather support the same critical behavior 
as for the pure system. 

In the S =1/2 case the data are not fully consistent 
with the expected RS exponent rf^s = 2, shown as a 
dashed line, but are characterized by an increasing slope 
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FIG. 3: Finite size analysis of the inverse dynamical exponent 
1/z', as function of the inverse chain length 1/L for S — 3/2 
at weak disorder 8 = 0.3. In the RS phase one would expect 
1/z' -> 0. 



as L increases in the log-log scale, which is an indication 
that the asymptotic regime has not yet been reached. 
This slope in any case is clearly larger than 1. Notice 
that for the calculation of the bulk exponent we use a 
combination [see Eq. of correlation functions be- 

tween an edge spin and a spin in the middle of the chain. 
As the distance between these spins is L/2, it is plausi- 
ble that Cl is plagued by stronger finite size corrections 
compared to the surface-surface correlation function C S L . 



B. Dynamical exponent 

We consider next the dynamical singularities in the 
system, which are related to the integrated probability 
distribution, 0(e), of the smallest energy gap e. This 
quantity, for small e, behaves as 



O(e) 



-1/z' 



(5) 



If z' < 1 quantum fluctuations dominate and the true 
dynamical exponent is z = 1, as in the pure system. If 
z' > 1 the disorder and quantum effects compete and 
the singular properties of the system are governed by a 
conventional random fixed point with z = z' . In the RS 
phase z = z' — ► oo. 

Figure [21 shows plots of In O(e) versus lne for 5 = 
3/2 (a) and S — 1/2 (b). The different scaling behavior 
for the two values of the spin can be seen by comparing 
the data for the same strength of disorder (5 — 0.3). 
In the 5 = 1/2 chain the distribution becomes broader 
by increasing the chain length, while in the 5 = 3/2 
chain In O(e) tends to have a finite non- vanishing slope 
for small e, implying that 1/z' > 0, a conclusion which 
is at odds with a RS phase. 

For sufficiently strong disorder (5 = 3, however, also 



the 5 = 3/2 spin chains shows a broadening gap dis- 
tribution (Fig. |2Ja)). The inset of [3 presents a graph 
of In O(e) plotted as a function of the rescaled variable 
lne/L^, with ip — 1/2, as expected in the RS phase 
[l6[ . Similar results were also found at 5 = 2. At such 
strong values of disorder, due to difficulties in the conver- 
gence of the DMRG method, we restricted ourselves to 
exact diagonalization data of rather short chains L < 10. 
The asymptotic behavior is however reached already for 
rather small chains, as the good collapse of the data of 
the inset of Fig. [21 demonstrates. 

We have also calculated the dynamical exponent from 
the distribution of the end-to-end correlation functions 
C£, which is expected to scale as 



p e {c s L ) ~ (ci) 



1/z' 



(6) 



Figure |21 shows a plot of the exponent l/z'(L) vs. 1/L, 
for 5 = 0.3 obtained from the analysis of the gap and of 
the correlation function (Eqs. JSJ and ©). Both esti- 
mates lead to a finite value z' * 0.13. One may wonder 
whether this is a genuine effect or just due to a very slow 
convergence to a possible asymptotic behavior consistent 
with a RS phase 1/ z'(L) — > 0. On the view of our numer- 
ical results we believe that the latter scenario is rather 
unlikely, as it would require also a non-monotonic scaling 
in L of l/z'(L) for the data obtained from the gap. 

In a similar way we estimated the dynamical exponent 
for other strengths of disorder. For 5 = 0.8 we find z' = 
0.8(2). At 5 — 2 and even stronger disorder the scaling is 
consistent with a RS behavior. The borderline between 
the RS phase and the weak disorder phase seems to be 
located at around 5=1. 



C. About the Harris Criterion 

The numerical results for the correlation functions and 
for the gap distribution strongly suggest that at weak dis- 
order the 5 = 3/2 chain has the same critical properties 
as in the pure case, i.e. the disorder is irrelevant |18| . 
For sufficiently strong disorder the DMRG results show 
also that the system is in a RS phase, therefore there is a 
transition induced by the strength of disorder. The anal- 
ysis of both the gap and correlation functions indicate 
that the region of irrelevance extends up to 5 = 1 . 

According to the Harris criterion [19j disorder is rele- 
vant if an appropriately defined correlation length [2£j in 
the pure system diverges with an exponent v such that 
2/u > d (here d = 1). In the present case, as the dis- 
order is coupled to the local energy operator, Sj ■ 5;+i, 
the correlation length £ entering in the Harris criterion 
is that associated to the strength of the dimerization, t, 
defined as the amplitude of the alternated modulation of 
the coupling constants: Jj = 1 + (— 1)H. In the limit 
t — » this correlation length diverges as £ ~ \t\~" ■ In 
the 5=1/2 case, v = 2/3 < 2 [2l|. implying relevance 
of weak disorder. 
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FIG. 4: Finite size scaling of the first (a) and second (b) gap 
in the dimerized 5 = 3/2 chain. Data collapse for the first 
gap is expected to work only in the disordered region (t > 0). 



In the 5 = 3/2 case there is no numerical result avail- 
able for the critical exponent, v, therefore we have ana- 
lyzed the scaling behavior of the smallest gaps as function 
of t. Scaling plots are given in Fig. 0]and as a compari- 
son the same quantities for the S = 1/2 chain arc shown 
in Fig. The gaps are expected to scale as: 



s k (L,t) 



(7) 



with Ek a scaling function and where k = 1,2 labels 
the first and second gaps [23|. As shown above for the 
5 = 3/2 at weak disorder and for the pure system the 
surface is ordered at t < 0, i.e. the correlation function 
between the edge spins is finite in the limit L — > oo (for 
L even). There are several known examples of models 
with spontaneous surface order [2^, as the one observed 
the 5 = 3/2 chain. In these cases one expects a localized 
low energy excitation at t = with a gap that vanishes 
faster than 1 /L because the ground state is degenerate in 
the thermodynamic limit, i.e. e\ ~ L~ Zs with z s > 1. In- 
deed this can be seen in Fig. 0]for the smallest gap in the 
5 = 3/2 case, e\. In absence of dimerization this quantity 
scales with an exponent z s s» 1.36 > 1, as shown by the 
scaling collapse at the point t = in Fig. |3{a). The fig- 
ure also shows that in the surface disordered region t > 
a reasonably good data collapse for e\ is obtained with 
the choice L 1 /" 3 /; with v s = 2. Notice that in the surface 
ordered region t < no scaling collapse can be observed. 
All the other gaps in the 5 — 3/2 are non-localized and 
scale with z = 1, as illustrated in Fig. 01b) for £2, and 
their scaling form involves the bulk correlation length, 
£ ~ \A~ V 1 with v pa 1. This estimate for v is different 
from the value v = 2/3 in the 5 = 1/2 chain. We expect 
that the origin of this difference is due to the existence 
of a dangerous irrelevant scaling variable [24[, which is 
generally observed at first-order surface transitions [2^ 
as in this case. 

Our numerical results, summarized above, can be com- 
pared with theoretical predictions in Ref . 2FJ| ■ In the res- 
onating valence bond picture the low energy excitations 



FIG. 5: As in Fig. Hfor S = 1/2. 



of the open 5 = 3/2 chain are expected to be described 
by effective 5=1/2 edge-spin degrees of freedom, which 
are very weakly coupled to the rest of an 5 = 1/2 chain. 
In the thermodynamic limit the edge spins are expected 
to be decoupled, the effective coupling being logarithmi- 
cally small, as 1/lni. In this picture the smallest gap 
(at t — 0) scales as e\ ~ 1/LlnL, whereas the other gaps 
scale in the same way as for the open 5 = 1/2 chain. 
Our numerical findings disagree with this prediction at 
two points. First, the localized gap has a faster size- 
dependence, our numerical results can not be described 
by a logarithmic correction. As a matter of fact for sizes 
we used in the calculation the effective exponent, z 8 , is 
monotonously increasing with the size and there is no 
sign of reversing this tendency. Data of previous numer- 
ical calculations in Fig. 5 of Ref. [27| show also strong 
deviations from a pure logarithmic correction. Our sec- 
ond disagreement concerns the second gap at t = 0, the 
value of which according to conformal invariance should 
be: 



£ 2 (-£<,0)L = TTV S X S 



(8) 



in the large L limit. Here v s — 3.87 is the sound ve- 
locity Q for the 5 = 3/2 chain and x s is the surface 
anomalous dimension, which for the 5=1/2 chain is|2S| 
x s = rj s /2 = 1. For the 5 = 3/2 chain the results in 
Fig. 01 are considerably smaller, our estimate is around 
x s = 1/2. Thus there are very probably new operators 
for the open 5 = 3/2 chain, which are not present in the 
5=1/2 chain. 

As a comparison we plot in Fig. [3] a similar scaling 
collapse analysis for the first two gaps in the 5=1/2 
case, using a similar range of system sizes as in Fig. 0] 
Notice that the best collapse is obtained with a scaling 
variable L t, implying v w 0.71, not far from the ex- 
pected value v = 2/3, a difference which can be imputed 
to logarithmic corrections. 

The analysis at finite t reveals that the dimerization 
couples differently to the 5 = 3/2 chain with respect to 
the 5=1/2 chain. We argue that the appropriate cor- 
relation length to be used in the Harris criterion [T^. |20| 
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FIG. 6: Schematic phase diagram of the spin S AFH chain 
as function of the strength of disorder 8. The RSs denotes a 
random singlet phase where the relevant degrees of freedom 
are effective 5-spins. The WD and ID are the region at weak 
and intermediate disorder. 

for the 5 = 3/2 chain is £ s that associated to the smallest 
gap and thus with v s = 2 the disorder is marginally irrel- 
evant. An accumulation of weak bonds along the chain 
may result indeed in an effective cut of the system. The 
resulting scenario is that of weakly interacting segments 
of finite length in which the relevant length scale is £ s , 
that associated to a chain with open boundary condi- 
tions. (We note that a similar scenario is used for the 
renormalization of random 5 — 1 chains, in which the 
effective coupling between spin- 1/2 degrees of freedom 
is exponentially small[llj.) This scenario therefore pro- 
vides a plausible explanation for the observed irrelevance 
of weak disorder for the random 5 = 3/2 chains. 

III. CONCLUSION 

We conclude by suggesting a general phase diagram 
for the spin-S AFH chains as a function of the strength 
of disorder <5. This phase diagram, shown in Fig. H3 is 
obtained by combining DMRG and SDRG results known 
for the 5 = 1/2, 5 = 1 and 5 — 3/2 cases with some 
general arguments for higher 5. First of all, differently 
from the 5 — 1/2 case, where any amount of disorder is 
known to lead to a RS phase, we expect that all chains for 
5 > 1/2 have a region of irrelevance, the weak disorder 
(WD) region, where the system is either critical (half- 
integer 5) or gapped (integer 5*) . This was already known 
for 5 integer. Our results for the 5 = 3/2 chain suggest 



that it may be true also for non-integer spin. 

At strong enough disorder all 5 are in a RS phase com- 
posed of effective 5 spins, which we indicate as RSs- For 
the S = 1 case there is only one such phase (RSi), while 
a recent SDRG study for 5 — 3/2 ^1 predicts two dis- 
tinct phases: RS 3 / 2 and RS^- We found no signatures 
of this spin reduction transition in our DMRG calcula- 
tions, which could be anyhow difficult to detect as the 
RS1/2 and RS3/2 have identical critical exponents. There 
is also the possibility to have, for higher spins (5 > 3/2), 
a sequence of multiple RS phases of different nature. Pos- 
sible phase diagram for higher 5 have been recently dis- 
cussed |13|. We have indicated these transition lines in 
the phase diagram of Fig. [ft] which, for the time being 
cannot be supported by numerical results. Finally, the 
WD region may be separated from the RS phase(s) by an 
intermediate disorder (ID) region where exponents vary 
continuously with 6, as observed in other models |25| . 
For the gapped (integer S) case this would correspond 
to a region of Griffiths singularities. The qualitative na- 
ture of the disorder induced cross-over phenomena for the 
5 = 1 chain is still debated (see in Refs. l29ll3(l ). We will 
discuss this issue in a future publication |31j . 
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Note added - After this paper was submitted we be- 
came aware of a recent work of Seguia et al. |32j |. who 
performed a SDRG analysis of the disordered 5 = 3/2 
AFH chain. Their calculation differs from that of Ref. 
|12| as the decimation scheme keeps into account more 
states, thus it is expected to be more accurate. While 
the results of Ref. [13] indicated the existence of two 
RS phases with different effective spins (SLff = 1/2 and 
5 c fF = 3/2), in the work of Seguia et al. [32j the renor- 
malization flow indicates that the weak disorder is an 
irrelevant perturbation of the system. This is in agree- 
ment with the conclusion of this paper. 
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